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Abstract 

A Slavnov- Taylor identity is derived for the gluon polarization tensor in hot QCD. 
We evaluate its implications for damping of gluonic modes in the plasma. Applying 
the identity to next to the leading order in hard-thermal-loop resummed perturba- 
tion theory, we derive the expected equality of damping rates for static transverse 
and longitudinal (soft) gluons. This is of interest also in view of deviating recent 
reports of 74 (p = 0) 7^ ^i{p = 0) based on a direct calculation of ji(p = 0). 
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1 Introduction 



Much interest is devoted to the physics of a the deconfinement phase of 
hadronic matter (quark-gluon plasma, QGP), with both strong experimental 
and theoretical research going on. A good theoretical understanding of this 
new phase of matter is urgently required, in particular, in order to identify 
worthwhile observables of this new phase of matter. Insight can be gained in 
the framework of 'Hot QCD' supplemented with a perturbative approach and 
important progress has been made during the last decade [1]. In particular 
thermal effects are known to alter soft modes in an important way, chang- 
ing dispersion relations, and allowing for Landau damping as well as for new 
modes of propagation for both fermionic (quarks) and bosonic (gluons) quasi- 
particles [2]. The collective dynamics involved can be consistently described 
within Hard-Thermal-Loop (HTL) resummed perturbation theory [3,4]. 

In this letter we are interested in damping of collective gluonic excitations 
as one of the most important characteristics of the plasma dynamics. While 
on-shell damping is known to be absent to HTL order g 2 T 2 a non zero damp- 
ing rate is expected to arise at order g(g 2 T 2 ). In early work by Braaten and 
Pisarski [5] a finite, gauge independent result for the damping rate of static 
transverse gluons 7i(c<j = m g ,p = 0) has been obtained, as one of the im- 
portant applications of the HTL-resummation scheme. Note that gauge in- 
variance also follows from a general theorem [6]. Also, in the static regime, 
Ji(mg,p = 0) = 7 t (m 9 ,p = 0) has been mentioned in [5], as is expected since 
longitudinal (1) and transverse (t) degrees of freedom can not be distinguished 
at zero spatial momentum. However, a complete understanding of the damp- 
ing behavior in the plasma is not yet available. In particular, for non-zero 
values of the spatial momentum singular results are obtained for transverse 
and also for longitudinal gluons [7], which originate from the absence of a 
chromo- magnetic mass in the HTL resummed gluon propagator. Moreover, 
also the analyticity of the damping rates around zero momentum has been 
questioned recently [8]: Expansion for small p of the imaginary part of the 
gluon polarization tensor has been reported to lead to 7 t (ra g ,0) ^ •ji(m g ,0) 
and even to a quadratic divergence in r yi(m g ,p = 0). 

In this letter we derive a Slavnov- Taylor (ST) identity for the gluon polariza- 
tion tensor in Coulomb gauge. After discussing some general consequences of 
it, we apply it to the next to leading order contribution to the gluon polariza- 
tion tensor TP", which might be helpful in order to further clarify the physics 
of damping for gluonic modes in the plasma. From it we easily obtain well 
known transversality properties as well as new constraints in particular con- 
cerning the imaginary part of TP" related to damping. The expected equality 
7i(m 9 ,0) = 7 t (m s ,0) can also be derived. We report on these findings in the 
next two sections and further implications are discussed in the conclusions 
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section. For convenience, the derivation of the ST identity is deferred to the 
Appendix. 



2 Slavnov- Taylor identity and its consequences 



Within the imaginary-time formalism of hot QCD, ST identities are derived 
from the BRST invariance of the Lagrangian in the same way as in vacuum 
theory. For a discussion of dynamical quantities continuation to Minkowski 
space should be done in the usual way. We briefly summarize the procedure 
in the Appendix. 

Most important to the dynamics of gluonic fields is the identity involving the 
gluon polarization tensor IP". In Coulomb gauge it reads: 



P V W V {P) = - 5^P 2 - P»P V + W U {P) W{P) 



with Hg(P) related to the Coulomb-ghost self-energy part Tl g (P) through 
p l U g (P) = n 9 (P).Q The derivation of the result (1) is explained in the Ap- 
pendix. Note that it is an exact identity, valid for arbitrary momenta P and 
to all order in perturbation theory. 



^From Eq. (1), other identities are immediately obtained that will be useful 
for the discussion below. First contracting also the remaining Lorentz index 
we obtain 



P fl W(P)P u = -P^W U (P)W Q (P) 

'IT" + 1 

Concerning the imaginary part of IP"(P), sinceQ ImH v {P) = 0, we obtain 



[p V - p»Pv\ KK + njn^n; . (2) 



P v IwR^{P) = -W g (P)ImU^{P) . (3) 
P^ImU^(P)P v = W g W g ImU, u . (4) 

4 Unless otherwise stated, throughout in this letter, Greek letters, fj,, u, take 
0,1,2,3 while Ratin letters, i, j, take 1,2,3. 

5 This can readily be seen in real-time formalism [4]. The building blocks of per- 
turbation theory, i.e., the propagators and vertices, and then also the self-energy 
part, are 2x2 matrices in "thermal space." ImW g {P) here is proportional to the 
(1, 2)-component of (IL g (P))ij = 1,2). Since, in Coulomb gauge, the ghost 
propagator-matrix (A g (P))^ = 1,2) is diagonal in thermal space, (JI g (P))ij is 
also diagonal and (U g (P))i2 = 0. 
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Note again that, as Eq. (1), also Eqs. (2) - (4), are general, valid for arbitrary 
momenta and to all orders in perturbation theory. 

Eqs. (1) - (4) contain important informations. At leading and next to leading 
order several of them have been established through direct calculation. Con- 
sider the soft momentum region, P M = (u, p) ~ gT; from Eqs. (1) and (2) we 
readily rederive two well known results from power counting arguments: 

(1) At lowest nontrivial order of effective (HTL-resummed) perturbation the- 
ory [3,4] W{P) (= 8W V (P)) = 0{g 2 T 2 ) so the left-hand side (LHS) of 
Eq. (1) is of 0(g 3 T 3 ). On the other hand W g (P) = 0(g 2 T) (no HTL- 
contribution in amplitudes with external ghosts [3,4]) and the right-hand 
side (RHS) of Eq. (1) is of 0(g A T 3 ). Therefore, to the HTL accuracy, 
P V W U {P) ~ P V 5W V {P) = 0. 

(2) Similarly, at leading order, the LHS of Eq. (2) is of 0(g A T 4 ). In the 
second line of Eq. (2), II£(P) = 0(g 2 T), W V (P) = 5W V (P) = 0(g 2 T 2 ), 
U g (P) = 0{g 3 T 2 ). Therefore, RHS/LHS = 0(g 2 ), and 

P v n^(P)P„ = (5) 

holds to relative order g 2 . A more precise statement is derived below, cf. 
Eq. (8). 

We now proceed exploiting the results Eqs. (1) - (4) in more detail. We focus 
on the soft momentum region and the next to leading order contribution to 
IP"(P) from which the leading damping behavior will be deduced. Following 
standard notations, we write W V (P) to this order 5W X {P) + TT' 1 (P) with SU 
the HTL-contribution and *I1 the relative order g correction. The diagrams 
that yield leading contributions to the imaginary part of *IT are shown in Fig. 1 
below. 

In Eqs. (1) - (4) we need an explicit expression for the ghost contribution. It 
is sufficient to calculate H%(P) to lowest nontrivial order, which we carry out 
in Appendix (cf. Eq. (32)): 



■ o 2 NT ■ 

W g {P)^-5^ 9 -^p\ (6) 



where p l = p l /p. It is worth mentioning in passing that U^(uj,p = 0) = 0. 
This is because, for p = 0, the summand/integrand of Eq. (30) in Appendix, 
which is valid for arbitrary (Euclidean) four- momentum Pe, is odd in the 
summation/integration variables Ke- 

First we use Eq. (6) in Eq. (1): To next to leading order under consideration 
we can write 5W V for W v on the RHS of Eq. (1) and using 5W U P U = [3,4] 
as well as (5) we find 
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^•n-(P) = ^^A-(f>)(i^p)" (7) 

with the inverse plasmon propagator [3,4], 



*AT\P) = p 2 - 5U 00 (P) = p 2 + 3m 



2p 



where m 2 = N(gT) 2 /9. Note that the close connection with the inverse prop- 
agator is of course apparent in Eq. (1). Also Eq. (7) is valid for any value of 
A, the Coulomb gauge parameter, Eq. (16) (in Appendix). We have confirmed 
Eq. (7) through direct calculation of the LHS. 



Incidentally, from Eq. (7), P M *rP(P)el r) (p) = 0, (r 
vector e\ r \p) perpendicular to p and arbitrary (soft) P M . 

Using Eq. (6) in Eq. (2) we find 



1,2) for any three- 



p/n^(P)p„ = 



'g 2 NTuo\ 



l-—p i 5n ij (P)p> + 0(g) 



(8) 



which is of 0(g 6 T 4 ) when P^ is soft, thus specifying somewhat the result under 
point (2) above. 

Finally, making progress towards the discussion of damping in the next sub- 
section, we note that, from Eq. (7), the imaginary part of P^TP^P) is related 
to the Landau damping contribution in the HTL-propagator. We have: 

Im \P;W»{P)} = ^9 2 NTm 2 g Q ' 9(p 2 - u 2 ) , (9) 

and therefore for time-like momenta, 

Im [P/IP^P)] = (P 2 > 0) . (10) 



3 Damping in a gluonic Medium 



Soft gluon damping rates relate to the imaginary part of the gluon polarization 
tensor U^, leading damping arises from the relative order g correction TP". 
The relevant diagrams are shown in Fig. 1 where however the Coulomb-ghost 
loop does not contribute to I mil 111 '. HTL resummed propagators and vertices 
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have to be used since the dominant contribution arises from configurations 
with soft loop momentum. For the following presentation we restrict to strict 
Coulomb gauge (A = in Eq. (17) in Appendix). 

For gluon damping only the transverse and longitudinal components of TP" 



*n t (P) = ^(p)*rF(p) , ^'(p) = s» - py , 

*U l {P) = -*U 00 {P). (11) 

will be relevant. In the vicinity of u = u> t /i(p), the mass-shell for quasiparticles 
to HTL-accuracy, the respective inverse propagator component can be written: 



G t/i - TZ t/ l(p) (w - ut/i) + iIm*U t/ i(u = u) t/h p) , 

x 3m^ f 2 - (cj - p 2 ) 2 p ^o 
Zt iP)= ^t-f) ' 



Zf\p) = 3p 2 



w 2 



uji(ujf-p 2 ) 3ui 



+ o 2p 2 



m g 



lYmm. this, the damping rates for longitudinal and transverse gluonic modes 
are obtained as: 



ltMp),p) = -Z t (p)Im*Il t (uJt(p),p) --^—Im*U t (m g ,0), 



2m g 



7 ( w , (p) ,p) = +Zi (p) im *n, ( w , (p) , p) ^ im *n, K , o) , (12) 



with m g = Wt/iip — 0) the plasma frequency. 

The direct calculation of 7^, for general p, has not yet been accomplished, 
however a leading logarithmic singularity has been extracted for non-zero mo- 
mentum in both 7 t and 72 [7]. At zero momentum the calculation of % is well 
documented in the literature but note that calculating 7; requires to expand 
*Hi(P) to order 0(p 2 ) - a formidable and subtle task. A recent controversy 
[8] suggests that approximations should be handled with care and observing 
general constraints. 

We are now in a position to show that our results so far leads to 7;(m 9 , 0) = 
7t(m 9 , 0) and thus a consistent calculation necessarily provides it. On the gluon 
mass-shell, P 2 = uo 2 — p 2 > 0, we obtain from Eq. (5) and Eq. (10) 

Im*U 00 (P) = -Im*Ui{P) = {^j 2 p l fPIm*U ij {uj, p) . (13) 
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Specializing to p = and using the fact *IT J (u;, p = 0) oc 5 tj we find 

fflm^iu, 0) = \ Im*U u (uj, 0) = Im*U t (uj, 0) . (14) 



Finally, using u>i/t(p = 0) = m g and Eq. (13) with (14) in Eq. (12), we derive 
that 



holds, as is expected on physical grounds. 

Finally it is interesting to observe what consequences arise if we assume, on 
physical grounds, that 7;(m 9 , 0) = 7 ff (m g , 0) to all orders in perturbation the- 
ory. From Eqs. (11), (12) this assumption leads to 



4 Conclusions and remarks 



The Slavnov- Taylor identity for gluon polarization tensor we report on in 
this letter, Eq. (1), allows to obtain various constrains. Applying it to the 
next to the leading order in HTL-resummed perturbation theory, the identity 
Eq. (7) results, which is relevant to leading-order damping rate of gluonic 
modes. As has been discussed, for zero three momentum the expected equality 
Jt{ m g, 0) = ■Ji(m g , 0) can be derived with the help of the identity Eq. (7). One 
of the advantages of the above procedure is that the explicit expansion of 
*U 00 (u}i,p) around p ~ 0, which has been found to be troublesome, can be 
avoided in this derivation. Moreover, concerning the direct calculation of 7z, 
the present work indicates that dealing with singularities in intermediate steps 
of the calculation and necessary changes in the integration variables should 
be made in a way consistent with the identity Eq. (7). 



li{m g ,0) = 7t(m p ,0) 



(15) 




and using this in Eq. (4) one finds 



lim \n» [IT; + 25 ut p { ] ImU^ v = . 
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Appendix Derivation of (1) and (6) 



We start in the imaginary-time formalism continuing to real energies at the 
final stage. The QCD action in Coulomb gauge reads 



S= J dx Jd 3 xC(x) = jd A xC{x)\ , 

V T / 



= -\Fr(x)Fr(x) - ^(V • A a (x)f + S^U^D^AMx) , 
F? = d»A» a - d»A» + gC abc AZA» . (16) 

Here D^ ab (A) = 5 a bd/dx^ — gC abc A^(x) and fj a and r\ a are the Coulomb-ghost 
fields. For notations related to the color space, we follow [9]. 

The bare gluon propagator A^ u (Pe) reads 

A^(P E ) = 5^5 U3 J F + <PV°4 + A^4^ , (17) 

p2 p 2 p 4 



where P£ = (p , p) with p = 27rTn (n — ± 1, ±2, ...), and the bare ghost 
propagator reads A 9 (p) = 1/p 2 . The form of gluon-ghost vertex can be read 
off from C 3 gC abc (d J fj a )r] b A J c , 

V g = -igC ahc j? , (18) 



where p is the outgoing momentum carried by i] a . 
The generating functional reads 

Z[J,H] = J VAZVrj a Vfj a exp S + J d 4 x{J^x)A^x) 

L T 

+^a(z)&(aO}] , (19) 

where, the functional integral is to be performed with periodic boundary con- 
ditions for all fields (including ghosts [10]), A£(x = 0,x) = A£(x = 1/T, x), 
etc. 

The action S (or C), Eq. (16), is invariant under the BRST transformation 
[9]: 

5A» = D^ b (A)r] b 5( , 5r] a = ~ V • A a 5( , 5 Va = 9 -C ahc r) h r] c 5( , 
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where 5( is a Grassmann-number parameter. Using this fact in (19), we obtain, 



Jd 4 zB(z)Z[J^,^] = 0, 



B(z) = JZ(z)D» ab 



S\ 8 1+ , s d 5 



&(z) A sav 'd*5J*(z) 
+ 2 UabcU(Z) 5£ b (z) 5Uz)- 



(20) 



Computing 



5 5 
5X(x) 5Uv) 



J d 4 zB(z) InZ[J,£,£] 



j=?=C=o 



by using (20), we obtain 



1 d 



5 ab d%G g (x -y)- gC acd (A^)Vc{x)n b {y)) + ^-—^{x " y) = ° ' (21) 



Note that the third term of Eq. (20) does not contribute to (21). The ghost 
propagator G g and the gluon propagator G pu in Eq. (21) are defined, respec- 
tively, through 



SabG g (x -y) = (r]a(x)r]b(y)) = 
5 ab G» v {x-y) = {A»{x)Al{y)) 



5\nZ 



5\nZ 



j = | = « = o ' 



5J£(x)5Jtty) 

(Aj(x)r) c (x)fj b (y)) is the gluon-ghost three-point function. 
Now, we call for the Schwinger- Dyson equations: 



&»(x-y) = A^(x-y) 

- J d 4 ud 4 v G" p {x - u)U pa (u - v)A au (v - y) , 



SabG g (x -y) = S ab A g (x - y) 
+gC ade d 4 z 



dA g (x - z) 
dz i 



{Ai{z)rj d {z)fj b {y)} , 



(22) 



(23) 
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where U pcr is the gluon polarization tensor. It can readily be shown that, in 
the limit g — > 0, (21) holds, as it should be. Taking this fact into account, we 
substitute (22) and (23) into (21) and use the form (17) for A at to obtain 



5 ab J d 4 zU^(x-z)d^A g (z-y) 

T 

= -gC acd J d A i J dhiG-'Y^x - 

T T 

>^fl|A fl (e - z) + 5 4 (£ - z)V*\ {A^z)r, e {z)fkiv)) ■ (24) 



x 



The inverse propagator (G' 1 )^ is related to through (G" 1 )^ = (A^ 1 )^ 
IF* 1 . The Fourier transform of {A p d (z)r] c (z)f] b (y)) may be written as 



gCa^A^z)^)^))]^ = -t5 ab W(P E )G g (P E ), 



(25) 



where (W g (P E )) ab is related to a ghost self-energy part Tl g through 

p%{P E ) = MPe) ■ 



(26) 



Putting altogether in (24) and transforming to Fourier space, we finally obtain 



n^(p E )p% = -p 2 G g {p E ) [{s^p 2 - p»p£ + w»{p E )}w g {p E ) 



-^U^(P E )U g (P E ) 



(27) 



where use has been made of (26). Noting that G g (P E ) = [p 2 + U g (P E )] 1 and 
using W u = ft"' 4 , Eq. (27) may be 'solved' as 



P»W»(P E ) = - 5^P E -P£P E + W»(P E ) W(P E 



(28) 



/,From this we obtain Eq. (1) in the main text after continuation to real ener- 
gies ipo — > uj + i0 + according to 

P" = (u>,p), P 2 e ^-P 2 =P 2 -uj\ 

rp^Po, p) -> HA° +a "°rp>, P ) , w g (t Po , P ) - HA°if>, P ) . 

(29) 



Finally we compute U.^(P) with soft P to lowest non-trivial order — one- loop 
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contribution — as required for the discussion in the main text. Using Eqs. (17) 
and (18), we have 



,(30) 



-9 2 NTJ2 

k 



(2tt) 3 (p - k) 2 
d 3 k p 1 - k i 



where f = (0,p), K E = (k Q ,k) and 5 ± ij (k) is as in (11). In Eq. (30), A(K E ) = 
1/K E for hard p, while, for soft p, the HTL-resummed *A(Ke) is substituted 
for A(Ke)- After carrying out the renormalization, one can easily see that the 
leading contribution comes from the first term (in the square brackets) with 
soft k region: 



fc o so ftfc V ' yv > 
1/T 

x J dre ik ° T Jd(p t ((,k)[l + n(0]e- (T , (31) 
o 

where pt(C, k) is the spectral function of the HTL-resummed transverse gluon 
propagator and n(() = l/(e^/ T — 1). On summing over ko and using n(() —T/( 
and / d(p t ((, k)/( = l/k 2 , we obtain 



IL$(P E )~-g i NT6* J 



d 3 k p l — (p ■ k)k % 



(2tt) 3 P(p-k) 2 

soft k 
. a 2 NT ■ 

-"^Se"^' (32) 

which is independent of A. Noting that in (32) has only spatial components 
and is independent of p , we have II£(P) = %(-Pb) (cf. Eq. (29)). 
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Fig. 1. Next to leading order contributions to the imaginary part of the ghion 
polarization tensor. All loop-momenta are soft and HTL-resummed propagators 
and vertices have to be used. 
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